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The Elliptic Inequalities in the Lunar Theory. 

By Ernest W. Brown, Fellow of Christ's College, Cambridge. 



[Continued from page 263.] 

The coefficients depending on the first power of the eccentricity having 
been obtained with sufficient accuracy, it is possible to proceed and obtain those 
depending on the square and higher powers of the same quantity. It is to be 
remembered that the value of c, the ratio of the two periods, was only the first 
term of an infinite series arranged in powers of the square of the eccentricity. 
Hence in obtaining the coefficients depending on higher powers of the eccentricity, 
it may be found necessary to have more terms of this series calculated. 

V. 

In the system of equations (4) put p = , there results 

2j2 9 U/, *. °. 2) A J, Ai-i, 2 + (»» °) A Ai-j-h-v + P. °] A J, i A -i-J-h-i\ = ° 

for all integral values of i, except i = 0. Now the suffix ± q denotes that the 
term to which it is attached involves the eccentricity raised to the positive power 
q, hence in this equation only even powers of the eccentricity will occur. 
Omitting terms which are of an order higher than Tq, we have, in accordance 
with the notation previously used, 

Ai, 1 = e i ' -4», -1 = £ i • 

We now put 

A = «,• + $aj 

instead of a jt where ha } denotes the new part depending on the eccentricity; 
aj has the numerical value before assigned and §a,j is evidently of the order Y\. 
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Hence q will only take the values , =b 1 . Performing these substitutions we 
obtain 

2AU, i, 0, 0)(o, + $a,)(a,_, + &*-«) + (/, i, 0, 1) w _ f + (i, », 0, - 1) «>/_,[ 
+ (», 0)S,{(% + ^a,-)(«i-j-i + &*#-i-i) + 2e A'-j-if 
+ [*, 0]2,{(«* + «a»)(a_i-i-i + Sa-i-j-i) + 2e/-i-i-i} = 0. 

If we omit the terms depending on Y\ the equation becomes the same as that 
used by Dr. Hill to obtain a t . 

It is necessary to see in what way c is involved in the coefficients (j, i, 
0,0), etc. Referring back to the general equations in Section I it was seen that 
c always occurred in one of the combinations 2i + cp, 2/ + cq. The coefficients 
(j, i, 0, 0), (i, 0) and [*, 0] will therefore be independent of c, while the two 
coefficients (j, i, 0, ±1) will only involve this letter in their numerators. As 
the last series of coefficients only multiply quantities which are themselves of the 
order Y%, it is sufficient for the purpose of obtaining a 4 + Sa t to this order to use 
the known part of the value of c which is independent of the eccentricity. 

As the quantities a % retain their previous signification, we obtain for the 
determination of ha t the set of equations 

2j|(i, i, 0, 0)(a;8aj_i + %_^a,.) + (j, i, 0, l)%-_« + (j, i, 0,— l)e}ej- t 

+ 2 (», 0)(a i _j_ 1 5o, + e/Ljf-i) + 2 [*, 0] (a.^^^a,- + e/L^^i)^ = 

omitting as before that for i = . 

The labour of calculation is conducted in a way similar to that which was 
followed in the case of e t , s' t . The amount seems greater on account of the larger 
number of terms to be calculated ; it is, however, not necessary to push the results 
to so many places of decimals. In addition, the coefficients (/, i, 0, 0), (i, 0) 
and [i, 0] are given by Dr. Hill in Vol. I, pp. 245-6 ; his notations for them are 
respectively (j, i), (i) and [*] . The terms in e i( s' t can be calculated exactly at 
the first approximation by using the values previously given for these quantities, 
and in obtaining the second approximation it is necessary only to find the 
changes in 5a 4 . I recall that a f = 1 when i = 0, and therefore for that value of 
i, Sai = 0. 

The values of 8a t are given below to eight places of decimals. The calcula- 
tions were made to nine places, and to that order a third approximation was 
found to be unnecessary. 
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tax = + .03938 170 Yl, Sa_ 1 = +.01376 519 F 2 ,, 

8a 2 = + .00046 113 F 2 , 5a_ 2 = + . 00002 216 F 2 ,, 

So, = + .00000 473 F 2 , Sa_ 3 = + .00000 026 F 2 ,, 

§a 4 = +.00000 005 F 2 , 3a_ 4 = + .00000 000 F 2 . 

For the sake of completeness the value of 8a , the part of a which depends 
on the square of the eccentricity is obtained. We have (Vol. I, p. 131) 



xu 



(usf 
Substitute in this equation 



= (Z> 2 + 2«iZ>)m + |wj 2 (m + s). 



u = (a + Sa^^a, + Sat + e£ + eg-<}? i+ \ 



s={a + 5a )2i|«-i-i + ««-i-i + «-«-i^ + ^-i-i^l^ 1 , 

where a has the numerical relation to {(i/n 2 ) h given by Dr. Hill. Since the result 
of such a substitution ought to be an identity for all powers of £, consider the set 
of terms independent of % ±c . Those on the right-hand side to the order Yl are, 
after division by a + 8a , 

2 i [(a i + KH(2* + 1 + ™) s + W\ + fm 2 (a._*_ 1 + Sa^O^ 1 ]. 
Put f = 1 and let 

#= 2a] (2» + 1 + mf + 2m 2 } , 
SH= 2M\(2i + 1 + mf + 2m 3 } . 

The right-hand side then becomes Z7 + SH when £ = 1 . 

In the expressions for u, s just given, put t = t , i. e. £ = 1 except where it 
occurs to the index ± c : 

w = («o + 5a )2<{a< + <K + £*£" + e<^ _c h 
* = Oo + 5a )2i{«i + *«* + e«T + e«£~ e f , 

and substitute in *«/(«<> + 5a )(^) J . Taking out those terms as before inde- 
pendent of £ ±c , we obtain as far as the order F 2 , 

(a„ + Sa„) a ^ = (2a, + 23a, + 2 e <£ c + ^TtH^ + 2^ + 2^° + 2 £ ,£-«)- 3 

_ l / j_2jk _L26i_ /(! _ JL ?!£/-« . JL I?l^!l\ 

~(2^r s V 2 Sa 4 2 2a, ^ 2 2a>. "*" 4_ (2a,)V 

/ 3 2K 3 2*' 3 % r - c , 15M\ 

X \ 1_ ^~ 2a7~"2"2a:^ ~~ 2 2a> + 4 (2a,) 2 / 
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which gives 

1 /. 2SK 3 (2 £t - + 2eQ 2 + 42^.2^ 
(XatfV 2% "*" 4 (2a*) 3 " / 

for the terms independent of £ ±c . 

Let this expression be put equal to J + hJ where J= l/(2as) 3 .- Taking 
account of the value of x which is /« (1 + rrif\n % , we obtain 



Dr. Hill has obtained (Vol. I, p. 144) 

-=(£)'D^]' 

= +.99909 31420 (-^-) . 



Hence 

^o_J_ g J 1 3-ff 
a ~~ 3 * J" 3 * g 

A ?^b 4. (Se i + 2e t Q a + 42^-2^ _ 2fei|(2^+l+wi) 8 +2m^ 
~ 3 2a* 4(2a s ) 2 ~" 32% j(2* +1+ m) 2 + 2m 2 J 

Reducing to numbers by means of the known results, 

— °= — . 13311 28 IJ, 
a 

and therefore the common factor for the series in u, s to the order Y$ is 
a + fe = (+ .99909 31420 — .13299 21 Yjfif-^A . 

VI. 

The next step forward requires us to obtain the series of coefficients (i , ± 2) . 
In equation (4) put p = ± 2 successively. It will be seen that the terms are all 
of the order F 3 at least, and that to obtain them to this order we only need to 
know Si, si and the values of a t \ those of Sa t are not required. In this case q 
only takes the values 0, ±1, ±2, terms which involve a higher second suffix 
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than ±2 being omitted. It is convenient to denote A i>% and A tt _ 2 by f { and// 
respectively on the same plan as before. The equations for the determination 
of/ { and fL f become 

Xj\u, i, 2, 2)o,_^+(y, i, 2, o)o^'_ ( + (y, », 2, i) e>8 ;_ f 

+ 2 (t, 2)(a,_,_ 1 / y + «,_ i _ 1 e,) + 2 [», 2](a_ i _,-_ 1 .# + elc-,-^)} = 0, 
2,{(/- », -*, ~2, 0)a j _ i f j + (j-i, -i, -2, —2)a,fJ- t 

+ (/—•*! — *i — 2 » — !)% e i-i 
+ 2[-», -2](a,_ i _ l ^+e l _ i _ 1 a,) + 2(--t I -2)(a_ l _,_ 1 //+«L l _ J l l8 ;)f = 0. 

The first line of the second equation has been transformed by putting j — i for/ 
in order to make it uniform for calculation with the first equation. The two 
differ now only in their coefficients and can be dealt with together. The coeffi- 
cients in the first line of each equation are transformed in the same way as in 
IV. It is evident that the known value of c is sufficient, all the terms being 
of the order TJ. The terms depending on e 4 , s' t are fully calculated at the first 
approximation, and we proceed by methods like those used in obtaining the 
numerical values of 5aj . 

The results are as follows : 

A = + .09402 355 A- =+ .03180 170 

f f 

<=i = — .06517 276 -w= -r. 01564 642 

J-o -*o 

■^r =+ .00132 915 -&- = + .00428 597 
-to -«-o 

f f 

<=f = + .00000 174 -ifa = + .00004 843 
-to -to 

Ai = + .00000 003 -A- = + .00000 049 
•to -*o 

-^ = +.00112 370 "%r = + -00006 457 

•to -*o 

-A = + .00001 161 •%! = + .0000O 066 

f f 1 

-^P = + .00000 011 w = + .ooooo 001 

-to *o 
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The maximum number of approximations necessary to obtain these results 
was four. The common divisor of f_ x and/j' is 

(2 — 2c) 2 ] 2(2 — 2c) 2 — 2 — Am + m % \ . 

Asc=1.07...., the factor (2 — 2c) 2 causes this divisor to be small and /_ 1( f{ 
correspondingly large. As/_! and/j' are comparable in size with / and f' , the 
number of approximations necessary for a given accuracy is much increased. 
The extra labor involved is not, however, very great, as very few of the terms in 
the expressions are affected. 

When this divisor is expanded in powers of m. , a series whose early terms 
(in the case of our moon) gradually increase, results. For 

3 « 201 3 

c — l = m — -nr — nr — .... 

4 32 

..^_-i s (l + |m + ^ 7 ™« + ....) 



( c _ if ,„* 

accounting in some measure for the slow convergence of Delaunay's coefficient of 
the term with argument 2D — 21 (see Section IX). 

VII. 

It is now necessary to determine the parts of A it±1 which depend on the 
cube of the eccentricity : there are none dependent on the square of that quantity, 
as has been pointed out. The chief difficulty here arises from the fact that it is 
necessary to obtain a further approximation to the value of c. This latter is, 
however, the principal object in view in this section, and it was necessary to push 
the computations thus far in order to achieve it. 

In equation (4) put p = 1 and let q take in addition such additional values 
that the terms depending on Y\ may be included. Reference is made to Sections 
III and IV . Using the previous notation as far as it is available, the equation 
for p = 1 is 

Pj, i-Aj-i, 0-4;, i + Qj + i, «-4; + *,-o-4/, -l H" > *i 1 » 2)-a.j, a4/-*, i 

+ (/, », 1, — l)4.-^-i,-i+(i, l)(A-j-i,<Ai+A^-i-i,-i) 

+ [»» i] (^-i-i-i, A -i + A i-4-i-i-i. -i) = °- 

We now put as before 
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and now 

A, i = e« + $e { , A it _j = ei + 8e[ , 

so that 6i, s'i denote the values found in the first part of this paper of the order 
3^ and 8s it 8s'i denote the new parts to be found of the order Y%. Also let 8c 
denote the new part of the ratio of the two periods of the order Yl, so that 
c -f 8c represents its full value to this order. Let us for a moment denote by 
F t , the divisor 

(2i + pc + p8c) 1 2(2i + pc + p8cf — 2 — 4m + m 2 } 

which is common to all the coefficients. The coefficients (i, p), [i,p] have in 
addition the divisor (2* + pc + p8c). Multiply the equation by F t ; it will then 
be free from fractions involving c + 8c in the denominators except as regards the 
coefficients just mentioned. The equation has been already solved when 8c, 8s t , 
8si, 8ai and terms of higher order are put equal to zero. Hence taking the 
variations of c, s { , e{, a { and including such terms of higher orders as are required, 
we have, after division by F t (in which expression 8c can be put zero), for j>= 1 

{Pj^aj-i + 2(t, 1) «!_,_! fJa, + { Qj+i, fy + i + 2 [i, 1] a_ i _ j _ 1 }8sj 
+ {Pj, Mj^i + 2 (t, 1) 8a i _ J _ 1 \e j + { Qj+i, i8a j+i + 2 [*, 1] &»_,_,_! }e( 



+ [«,_, I (P A ^J) + 2a,_ i _ 1 | {(.", 1) Jty ] 3" 



+ (/, », 1, 2)/^_, + (y, », 1,-1) «(#_,+ 2 (., l) i ^ l / _ i _ 1 +2[t > 'l]e,/l < _ i _ 1 = 0. 

Here 2, is to be understood to stand before each term. There will be a similar 
equation for p = — 1 , 8s t and <$e'_« being a pair which are calculated together. 

The method of proceeding is by continued approximation as in the deter- 
mination of e t , si- We obtain 8ei, 84 for a sufficient number of values of 'i, 
excluding i = 0, in terms of 8e , 8s' , T 8c, Yl, the results as is evident being 
linear functions of these quantities with numerical coefficients. The results are 
then substituted in the two equations given by i= (or two equivalent equations) 
and we obtain two linear relations between the same four quantities. 

The number of terms to be computed is very large, but the following obser- 
43 
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vations will show that the amount of calculation is really smaller than at first 
sight appears. 

There are four sets of terms involving in combination (i) the letters a, he; 
(ii) the letters s, Sa; (iii) a, e, he; (iv) s,f. These can be dealt with sepa- 
rately at the first approximation. 

In set (i) we are going to determine the parts of he it hs'i which depend on 
he , &sl>. But the expressions are exactly the same as those given in IV for 
determining e t , e' t in terms of e 0t e' . The results are therefore to hand. 

In set (ii) the coefficients are the same as those in IV and are therefore 
known ; also as Sa u s u e' t are known, the set can be completely calculated at the 
first approximation in terms of F 3 . 

In set (iii) the new coefficients have to be found, a t , s t , e'i are known. 
These terms can also be completely obtained at the first approximation, ex- 
pressed in terms of T §c; as Sv comes out approximately .0027 F 3 , these terms 
are calculated to two fewer places of-decimals than those involving Y$. 

Set (iv) can be calculated completely at the first approximation in terms of 
F 3 , as all the quantities are known. In addition the coefficients (*, 1), [i, 1] 
have been found, while as 

(J, i, 1 , q) = (2/ + go) U + (2/ + qof M u 

and that L { , M t have been calculated, these coefficients are soon obtained. 

Having calculated as far as possible all the terms at the first approximation, 
we use the values of Set, 8si thus obtained for the second approximation. Set (i) 
will be the only one affected, and its coefficients have been found. These 
approximations will of course add to the coefficients of I^ 3 and Y Sc, but not to 
those of Se and Ssq, as is evident from the remarks on set (i). 

The coefficients s i} e' t have been expressed in terms of Y or s — s' . It now 
appears that Y has an additional part 3e — 8ei to be added to it. This really 
amounts to introducing a new arbitrary constant, or, in other words, to making 
£o — e' + Ss — 5eo the arbitrary constant instead of s — 4- This is inconvenient 
and causes unnecessary complications. Since we are at liberty to choose our 
arbitrary, it will be simplest to make it F , so that £F = 0. Hence the coeffi- 
cient of the principal elliptic term in the expression for the coordinate y will, in 
this theory, be 

a + Sa Y , 
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where a + 8a has the value which is found above expressed in terms of (/i/n*)* 
and Yq. We have then ^ = <^o • 

The values of <^, Sei found as above are given below, 

&•_!= + .02435 517 F? + .22567 789 $e + .76306 1 Y^c, 

Ssi — —.00877 697 I? —.08834 235 os —.24185 2 Y hc, 

$ El = + .02529 357 Y$ + .00215 229 Ss —.00006 1 Y Sc, 

Se'_i — + .00322 634 F 3 — .00128 260 o~e + .00079 8 Y Sc, 

&?_ 2 = + .00070 947 Y S + .00008 499 8e + .00003 7 F 3c, 

Sej, — +.00526 427 F 3 — .00049 147 &? — .00145 3 Y^e, 

Ss % = + .00048 131 TJ + .00000 620 Se — .00000 1 F Se, 

8sL z = + -00002 534 F 3 — .00000 076 &■„, 

&_, = + .00000 493 F 3 — .00000 019 Ss , 

Ses = + .00009 895 I? — .00000 333 Ss — .00001 Y 8c, 

$e 3 = + .00000 678 F<? —.00000 001 5e . 
Sei 3 = + .00000 031 F 3 , 
Se_ 4 = + .00000 006 F 3 , 
Sei = + .00000 137 F 3 . 

In order to obtain the numerical ratios of 8s , Sc to F,f, we may use the two 
equations of the system just considered for i= 0. As will be stated below, an 
equation of verification was, however, computed, — in this case, the second of 
equations (6) in the first part of this paper, with the necessary terms of order Fjj 
included. These three equations gave respectively 

+ .00425 226 F 3 + .02946 89 &■„ — 1.11898 7 Y Sc = 
+ .00902 093 Yi + .00129 30 Ss — 3.33860 6 Y hc = 
+ .18016 812 F> + 4.32005 66 $e + .987815 F Sc=0. 

The first of these combined with the third gave 

8c = + .00268 561 F<?, 

and the second with the third gave 2 instead of 1 for the last place of decimals. 

Whence 

&=„=&?/= -.04231 912 YI 
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Substituting in the table just given for the values of 8si, we obtain as the final 
values for these 

%^ = +.01685 40 -^-=—.00568 79 -fc- = + .0 25 20 23 -fej 1 = +.00328 28 

■H) -*0 -*0 -*0 

is V / fr § / 

%^= + .00070 60 4|- = +.00528 12 -^r = + .00048 10 '-^ = +.00002 54 

-M) *Q -*0 -*0 

%- 3 = + .00000 49 -^= + .00009 91 ^-= + .00000 68 -^ = +.00000 03 
-*0 -*o ■*<) ■'0 

M 

=4 =+.00000 14 

It is not my intention to carry the calculation of the coefficients to any 
further extent. Enough has been given to show how they can be obtained to 
any desired accuracy ; and the methods outlined above for making the computa- 
tions can be arranged so that the assistance of a computer of ordinary intelli- 
gence can be largely made use of. One of the chief advantages of these methods 
is the ease with which equations of verification can be employed. Before com- 
paring some of the results with the corresponding ones of Delaunay, I shall give 
a short account of the methods of verification used. 



VIII. 

Verification of the Results. 

There appear to be two principal sources of error in the class of computa- 
tions here treated apart from that produced by the use of logarithms. The first is 
the common one of making an actual numerical mistake of any kind whatever ; the 
second arises from a term omitted altogether. It is therefore necessary to have 
such equations of verification as shall give some security that errors proceeding 
from both these sources may be eliminated. Owing to the interdependence of 
the various approximations and, in the case of the higher coefficients, to the 
large number of terms to be computed, it is advisable to have at each step a test 
for errors arising from the first and chief source. For this object the equation 
most useful, on account of its simplicity, is the second of equations (3). It was 
used in all the calculations for coefficients depending on the square and cube of 
the eccentricity. In it the terms are arranged in the same way as those in the 
equations used for finding the results ; that is to say, the suffixes proceed in 
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the same order, and the coefficients depending on m, c alone differ. An error 
therefore of a sensible term omitted altogether in finding any pair of coefficients 
would probably also occur in the equation of verification. Such a case occurred 
once in these computations. For the coefficients depending on the first power 
of the eccentricity, one of equations (6) was used, as the computations of the 
coefficients of its various terms were simple, and as it served at one calculation 
to verify the four quantities of the same order, e it eLi, e[, e_i, while the second 
of equations (3) required two calculations, one for e 4 , el 4 and another for si, s_ { . 
It will be noticed, too, that the terms in equations (6) are arranged somewhat 
differently, a property which possesses some value in this connection. 

To detect an error of the second kind, reliance was placed on the computa- 
tion of an expression which has an entirely different form. In the equation 

9 m 2 
uDh + sD 2 w + Du.Ds — 2m (uDs — sDu) + -^- (u + s) 2 = G 

substitute the values of u, s, Du, etc., found and arrange the result in powers 
of F and £ ±c . The result will be of the form 

(^T^ =a+a '^ + c&^+^JW+ro+r*?^ +£~ 20 ) + • • • • 

Here a, a', /3, /?', y will be series in £ 24 with known numerical coefficients. As 
this equation should be true for all values of £, we should have 

° » = a + o'IJ, 0F O + /?'!?= 0, y = 



(a + 8a J 

to the degree of approximation used; moreover, since Y is arbitrary, in addi- 
tion we have 

(8=0, /?' = 0. 

Again, as these equations ought to hold for all values of £, we can put £ = 1 in 

each of them ; the verification depends on the results being identically zero. Its 

value arises from the fact that the expressions to be here calculated are all of 

such forms as 

2ia 4 , 2(2i+l)a u XiZi + l+cfsi, etc., 
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they are therefore totally different from any of those previously used. The 
higher coefficients as fo 3 , / 4 , etc., are, however, multiplied by large numbers in 
some cases, and therefore unless these are carried to a higher degree of approxi- 
mation than was before found necessary, this method fails to test quite fully the 
results obtained. For practical purposes, since the coefficients have been com- 
puted to a large number of places of decimals, the test appears to be sufficient. 

The Jacobian constant is obtained by putting^ = in the first of equations 
(3) . Let G now denote the initial value of this constant as found by Dr. Hill, and 
hO the new part depending on Y%. From the equation we obtain by means of the 
known values of Sa it 

— = + .34552 75 FJ. 
«o 

The value of SG obtained from the equation of verification above where £ is put 
equal to unity, differs only by one unit in the last place. This therefore gives 
the required verification for the series of coefficients Sa { . The coefficients f it /i 
are absent from these expressions. 

The case of these coefficients 8a t requires special mention. The two values 
of &G found from two entirely independent expressions differ only by one unit in 
the seventh figure. That is to say, they are the same within the limits of error. 
It seems very improbable that an error of one unit in the fourth significant 
figure of ha x or 8a_ lf or of the first significant figure in 8a 2 or &*_ 2 , should have 
been able to elude detection while running through two separate sets of equa- 
tions of verification. Yet this is roughly the amount of error required to make 
the results agree with those of Delaunay as shown in the next section. 

In order to verify the transformations from rectangular to polar coordinates, 
in which the method of ' special values ' was adopted, an extra ' special value ' 
was computed and the results were found to be accordant. 

All computations once made were gone through a second time. The average 
error made in the later portions of the work was about one in every four or five 
hundred figures. It did not seem to be confined to any particular class of 
operation. In using the new eight-figure tables of the French Government, 
extra care was exercised for the differences, and the chance of error thus 
diminished. 
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IX. 

Transformation to Polar Coordinates. 

Let v denote the difference between the true and mean longitudes. Gath- 
ering together the results of the previous sections and adding those given by 
Dr. Hill (Vol. I, p. 248), we have with Delaunay's notation for the arguments ; 
that is, 2D for that of the ' Variation ' and I for that of the principal elliptic 
term : 

r C ,° S ^ = 1 + (— .00718 00395 + .05314 689 FJ) cos 2D 
a + 8a ■ 0/ 

+ (+ .00000 60424 + .00048 328 FJ) cos AD 
+ (+ .00000 00325 + .00000 499 FJ) cos 6D 
+ (+ .00000 00002 + .00000 005 FJ) cos SD 

+ (—.49679 1802 F — .08463 82 FJ) cos I 

+ (—.09332 9284 F +. 01116 60 FJ) cos {2D — I) 

+ (+.00025 6338 F + .00598 71 FJ) cos (4D — Z) 

+ (+ .00000 2210 F + .00010 40 Fjj) cos (62) — I) 

+ (+ .00000 0025 F + .00000 14 FJ) cos (SD—l) 

+ (+.00134 2824 F + .02848 51 FJ) cos (2D + I) 

+ (+.00001 0769 F + . 00050 64 FJ) cos (4Z> + I) 

+ (+.00000 0079 F + . 00000 71 FJ) cos (6Z> + I) 

+ .12582 524 FJ cos 21 

— .04952 634 Y z cos (2D— 21) + .00118 827 FJ cos (2D +27) 

+ .00561 512 Y% cos (4D— 21) + .00001 227 FJ cos (42? +2?) 

+ .00005 017 FJ cos (6D—2Z) + . 00000 012 FJcos(6D+2Z) 

+ .00000 051 FJcos(8D — 21). 



rsm v 
a + <k 



+ (+.01021 14544 + .02561 651 FJ) sin 2D 

+ (+ .00000 57149 + .00043 897 FJ) sin 4D 

+ (+ .00000 00276 + .00000 447 FJ) sin 6D 

+ (+ .00000 00002 + .00000 005 FJ) sin 8D 
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+ (+ 1.00000 0000 F ) sin I 

+ (+ .20446 5776 F — -.02254 19 F 3 ,) sin (2Z> — 7) 

+ ( + .00036 0509 F + . 00457 52 F?) sin (4Z> — I) 

+ ( + .00000 1960 F + . 00009 41 F 3 ,) sin (6Z> — 7) 

+ (+ .00000 0023 F +. 00000 13 Yl) sin (SB — I) 

+ (+ .00159 6237 F + . 02191 96 F 3 ,) sin (2Z> + I) 

+ (+ .00000 9426 F + . 00045 57 Y 3 ) sin (4l> + I) 

+ (+ .00000 0069 Y + .00000 65 F 3 ,) sin (6Z> + 7) 

+ .06222 185 Y\ sin 2? 

+ .08081 918 Fgsin(2Z>— 27) + .00105 914 F* sin (2D + 27) 
+ .00295 682 F* sin (42) — 21) + .00001 094 Y\ sin (4Z> + 27) 
+ .00004 669 F| sin (QD — 27) + .00000 Oil F 3 sin (6D + 2?) 
+ .00000 046 F? sin ($D — 21). 

To find the longitude we can proceed as follows. The difference of the true 
and mean longitudes is denoted by 

v =. tan v — i tan 3 v + £ tan 5 v — .... 

Tan v is obtained from the above values and the terms are expanded in powers 
of F . To F° are attached the sines of the arguments 2iD, to FJ the sines of 
2iD ± I, and to Y\ those of 2iD, 2iD ± 21. The coefficients depending on FJ 
have already been obtained and discussed in the first part of this paper. As far 
as the order Yl let 

r cos v = So+tfi Y + S Z Y* 

rsmv=S!> + SlY + SlYl 

where S , /Si, etc., denote the terms involving the respective powers of F , etc., 
in the expressions just given. It is not difficult to prove that when these are 
substituted in the expansion for v, we obtain for the terms involving Yl in v 

S Si - SJ& , S S' (SI - Si 2 ) — Mi (go — W) 

si + s? ~*~ (si + wy "'" • 

The arguments ± 21 only appear in the numerators. The coefficients in this 
expression of cos 21, sin 21, and that in which 21 does not appear, are quickly 
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separated out and the results are series involving sines and cosines of the argu- 
ments 2iD. To each of these three sets of terms the method of special values 
for 2D can then be applied and the required coefficients of the terms in longi- 
tude obtained. 

The results for these terms are : 

+ r 2 [+ .12015 93 sin 2D + .30993 73 sin 2? +.08545 32 sin (2D— 21) 

+ .00209 70 sin 4,D + .00532 16 sin (2D + 21) +.01248 18 sin (4JD — 2?) 
+ .00002 70 sin QD + .00006 81 sin(4D + 21) +.00021 37 sin (62) — 2lj 
+ .00000 03 sin 8Z> + .00000 07 sin (6D + 21) +.00000 27 sin (SD — 2Z)]. 

The relation between F and Delaunay's e has been found in Section IV 
By means of it these coefficients expressed in-seconds of arc become 

+ 298". 959 sin 2D + 771". 132 sin 21 + 212". 610 sin(2Z>— 21) 

+ 5". 217 sin 4D + 13". 240 sin (2D + 21) + 31". 055 sin (4D — 21) 
+ 0".067 sin QD + 0".. 169 sin (4D + 21) + 0". 532 sin (6Z> — 21) 
+ 0". 001 sin SD + 0". 002 sin (6Z> + 21) + 0". 007 sin (SD — 21) . 

The part of the motion of the Lunar Perigee which depends on TJ nas been 
found above to be 

Sc= + .00268 561 Y$; 

and the relation between F and Delaunay's eccentricity, 

r = 2eX 1.00027 136. 
Also s / 1 ^»\ 1 s 

henCe t(±fo)= -.00994 29 e*. 

\ n at / 

The series given by Delaunay for this quantity is* 
675 „ . 31605 , , 1483665 



r 6 2 , 675 3 ■ 31605 4 , 



64 * ■ 512 1T 4096 



m? 



25291729 6 352038885 7 T 
+ 16384 mi ~ 1179648 "M"' 



*Comptes Rendus, Tom. LXXIV, p. 19. 
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where m x =.n' /n. In numbers the portion within the brackets becomes 
.00209 82 + .00441 42+. 00193 25 + .00084 82 

+ .00027 04 — .00000 39= .00955 96, 

which is considerably less than the previous value, the difference being greater 
than the last term but one of this series. Moreover, Delaunay's series seems to 
indicate that when estimation is made for the uncalculated portion, his value 
should be even less still.. The results, therefore, can scarcely be considered 
accordant. The difference amounts to 0".055 in the motion. 

The coefficients found above agree with what might have been expected 
from Delaunay's series, with one exception, the part of the coefficient of the 
variation which depends on e 2 . Expressed in seconds of arc, Delaunay's values 
for the terms e^m , e 2 m s , etc., in the coefficient are* 

217".976 + 59".839 +16". 331+ 3".852 + 0".777 + 0";137 — 0".049 = 298".863; 

and it would appear from this that the true value- was somewhat less, say 

298". 84 
the value obtained above is 298".96 



a difference of 0".12 

It should be pointed out that in Prof. Newcomb's comparison of Delaunay's and 
Hansen's Lunar Theories, f there seems no reason to expect that the former's 
value of the coefficient is erroneous by so much as one-tenth of a second. It 
may, however, have happened that there is an error in some other part of the 
long expression for this coefficient which may balance a possible error here, or, 
what is more likely, that the other terms if carried further would have given 
some considerable coefficient. The terms in e n converge very slowty up to the 
last given. 

It may be mentioned here, that in comparing the results for many of the 
coefficients with those given by Delaunay the differences are less than a hasty 
inspection would have caused us to estimate. If we examine the series arranged 
in powers of m , more carefully, and especially those which have a power of the 

*Mem. de l'Acad. des Sc, Tom. XXIX, p. 815. 
t Astron. Papers of Amer. Eph., Vol. I, p. 92. 
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eccentricity as a factor, it appears that in most cases the terms have a tendency 
sooner or later to change sign and the term before the change of sign has a com- 
paratively small numerical multiplier. Delaunay sometimes stops at this term, 
thus giving a false idea of the accuracy of the total result. 

The most serious difference is that in the value of the part of the Lunar 
Perigee just obtained, depending on e 2 and the ratio of the mean motions only, 
since, as Dr. Hill has pointed out,* the motion of the Perigee is capable of being 
found observationally with a high degree of accuracy. It is seen from the 
expansion reproduced above that Delaunay in this very slowly converging series 
has carried his terms.beyond the change of sign, and another change of sign in 
the uncalculated portion, which would be necessary to make the two results 
accordant, seems hardly likely to occur early enough to produce the difference. 
In either theory the difference, just mentioned, in the coefficient of sin ID would 
be sufficient to cause the discordance. As I have said in Section VIII, it seems 
hardly probable that a sufficient error has been made in the calculations under- 
taken to obtain Sa t . It might be surmised that the error arises from the set Ssi, 
Ss't The error required there would be several units in the third place of 
decimals, and again this seems improbable in view of the two sets of equations of 
verification. 

One other coefficient is worthy of special notice, namely that of sin (2D — 2?) . 
Delaunay gets for thisf 

130". 786 + 46".089 + 22".283 + 8".131 + 3".179 + 1".308 + 0".542 

= 212".318 

my coefficient being 2 1 2". 6 1 

so that the remainder of Delaunay 's series is + 0".29. The values of/ 4 ,/j' not 
depending on those of 5a 4 and vice versa, an error in one set does not involve one 
in the other. 

In a similar way may the parts depending on Y% be transformed to polar 
coordinates. The coefficient of sin I in longitude will be of the form 

aY + pY* 

* On the Part of the Motion of the Lunar Perigee, etc. Acta Math., Vol. VIII, p. 1. 
tMem. cit.,p. 823. 



338 Brown: The Elliptic Inequalities in the Lunar Theory. 

where a , (3 are numerical coefficients. In Delaunay's theory this coefficient is 

2e — Je 3 , 

and therefore the second approximation to T in terms of e will be given by 

<xF = 2e — ie 3 — ^f e 3 

(X 

This must be applied to all terms with arguments 2iD ± I when it is desired to 
express them in Delaunay's notation. 

Haveefoed College, June, 1893. 



